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The source of human infection with Middle East respi- 
ratory syndrome coronavirus remains unknown. Molecular 
investigation indicated that bats in Saudi Arabia are infected 
with several alphacoronaviruses and betacoronaviruses. 
Virus from 1 bat showed 100% nucleotide identity to virus 
from the human index case-patient. Bats might play a role 
in human infection. 


ince Middle East respiratory syndrome (MERS) was 

described in September 2012, over 90 cases have been 
reported worldwide, 70 from Saudi Arabia. The incidence 
of infection with the causative agent, a betacoronavirus 
(MERS CoV) (/), has not been determined; however, the 
mortality rate among those who received clinical care is 
265% (2). Although instances of human-to-human trans- 
mission have been documented between case-patients and 
others in close contact (including hospital patients sharing 
rooms, family members, and medical personnel), the sourc- 
es of infection for most patients remain unknown. Because 
of sequence similarities between B-CoVs identified in bats 
and those of MERS CoV isolated from humans, a bat reser- 
voir has been posited (3—5). Although neither detection of 
MERS CoV in bats nor contact of human MERS patients 
with bats have been reported, a role for bats in human in- 
fection cannot be excluded because contact can be indirect 
(mediated through another animal vector or fomites). 
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The Study 

In October 2012 and April 2013, three agencies col- 
lected samples from bats in regions where MERS cases had 
been identified (Figure 1). The agencies are the Ministry of 
Health of Saudi Arabia, the Center for Infection and Im- 
munity of Columbia University, and EcoHealth Alliance. 

During the October investigation, the team inter- 
viewed the family of an index case-patient in Bisha and 
collected samples from bats <12 km from his home, in an 
abandoned date palm orchard, and <1 km from his place 
of employment, a hardware store that fronted a garden and 
date palm orchard. Although neither family members nor 
employees recalled seeing bats, the team observed roosting 
bats and guano in abandoned wells and ruins within 12 km 
of his home and insectivorous bats at dusk in the garden 
behind his store. Over 3 weeks, 96 bats representing 7 spe- 
cies (Rhinopoma hardwickii, Rhinopoma microphyllum, 
Taphozous perforatus, Pipistrellus kuhlii, Eptesicus bot- 
tae, Eidolon helvum, and Rosettus aegyptiacus) were cap- 
tured in mist nets and harp traps, then released after visual 
speciation and collection of morphometric measurements; 
wing punch biopsy samples; blood; throat swab samples; 
and rectal swab samples or fecal pellets. Samples were col- 
lected into viral transport medium or lysis buffer. 

During the 3-week April investigation, fecal samples 
were collected on tarps laid out at bat roosting sites in and 
around Bisha, Unaizah, and Riyadh. Representative animals 
at each roosting site were captured, identified morphologi- 
cally, and released after wing punch biopsy samples were 
collected for speciation by DNA analysis. Samples were 
collected into cryovials. 


Bisha area 
Oct 2012 
Apr 2013 


Figure 1. Bat sampling sites and locations of home and workplace 
of index case-patient with Middle East respiratory syndrome, Bisha, 
Saudi Arabia. 
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DISPATCHES 


Table 1. PCRs and primers used in CoV detection* 


Nested fragment size, 


region (primer locations on Type of 
PCRs (reference) Primers, 5'3' the reference genome)t CoV (no.) 
Nested pan-CoV-| PLQ-F1, CGTTGGIACWAAYBTVCCWYTICARBTRGG ~400 nt, RdRp a-CoV (8), 
(6) PLQ-R1, GGTCATKATAGCRTCAVMASWWGCNACATG (18310-187450) B-CoV (1) 
PLQ-F2, GGCWCCWCCHGGNGARCAATT 
PLQ-R2, GGWAWCCCCAYTGYTGWAYRTC 
Nested pan-CoV-Il WT-COV-F1, GGTTGGGAYTAYCCHAARTGTGA ~430 nt, RdRp a-CoV (5), 
(7) WT-COV-R1, CCATCATCASWYRAATCATCATA (15260-15700) B-CoV (2) 
WT-COV-F2, GAYTAYCCHAARTGTGAYAGAGC 

WT-COV-F3, GAYTAYCCHAARTGTGAUMGWGC 
Hemi-nested RdRp- EMC-SeqRdRP-Rev, GCATWGCNCWGTCACACTTAGG ~230 nt, RdRp a-CoV (2), 
sequence assay (9) EMC-SeqRdRP-Fwd, TGCTATWAGTGCTAAGAATAGRGC. (15048-15290) B-CoV (1) 

EMC-SeqRdRP-Rnest, CACTTAGGRTARTCCCAWCCCA 

Hemi-nested N- EMC-SeqN-Fwd, CCTTCGGTACAGTGGAGCCA ~280 nt,N seq NA 
sequence assay (9) EMC-SeqN-Rev, GATGGGGTTGCCAAACACAAAC (29,549-29,860) 

EMC-SeqN-Fnest, TGACCCAAAGAATCCCAACTAC 
Nested Cll-pan- NM-CoV-2F1, ACWGTTCARGGICCWCCIGG ~355 nt, helicase B-CoV (2) 
CoV-III NM-CoV-2F2, GTTCARGGGCCWCCGGGNAC (17,060—17,410) 

NM-CoV-2R1, GGCAGCTGWGCWGGRTCICCNACRTA 
NM-CoV-2R2, AGCTGWGCWGGRTCGCCIACRTANAC 
Nested CIl-MERS- NM-HCOV-F1, GTGCTAAGAATAGAGCTCGCACT ~190 nt, RdRp B-CoV (1, 
RdRp NM-HCOV-F2, AGAGCTCGCACTGTTGCAGGC (15068-15249) MERS 
NM-HCOV-F2, AGAGCTCGCACTGTTGCAGGC CoV) 

NM-HCOV-R1, ACCCATAAGATGCGGATTATCAAC 

NM-HCOV-R2, TGCGGATTATCAACATCTTTGTAC 
Hemi-nested Cll- NM-NSeq-F-1, ACTTCCTTCGGTACAGTGGAGC ~170 nt, N seq NA 
MERS N sequence NM-NSeq-R-1, GGCACTGTTCACTTGCAATC (29545-29713) 

NM-NSeq-R-2, GGAGGTTCAGACATTTGGTCT 

upE and ORF 1b upE-Fwd: GCAACGCGCGATTCAGTT Upstream of E gene and NA 


real-time assays (8) 


upE-Prb: FAM-CTCTTCACATAATCGCCCCGAGCTCG-TAMRA 


ORF 1b 


upE-Rev: GCCTCTACACGGGACCCATA 
ORF 1b-Fwd: TTCGATGTTGAGGGTGCTCAT 
ORF 1b-Prb: FAM-CCCGTAATGCATGTGGCACCAATGT-TAMRA 
ORF 1b-Rev: TCACACCAGTTGAAAATCCTAATTG 


*CoV, coronavirus; MERS, Middle East respiratory syndrome; RdRp, RNA-dependent RNA polymerase; NA, not applicable; ORF, open reading frame. 
tPrimer locations are based on human B-CoV 2c EMC/2012, complete genome (GenBank accession no. JX869059). 


All samples were stored in liquid nitrogen and con- 
veyed to Riyadh for storage at —80°C before being trans- 
ported to Columbia University in New York in dry nitro- 
gen. The October 2012 shipment was inadvertently opened 
at customs in the United States and sat at room temperature 
for 48 hours before transfer to Columbia University; at ar- 
rival, all samples had thawed. The April 2013 samples ar- 
rived intact. 

Total nucleic acid was extracted from samples by us- 
ing the NucliSENS easyMAG system (bioMérieux, Dur- 
ham, NC, USA) and subjected to 8 PCRs with primers and 
protocols designed to amplify regions within the helicase, 
RNA-dependent RNA polymerase (RdRp), and nucleo- 
capsid or envelope proteins of CoVs (6-9). Products were 
sequenced and analyzed for similarity to GenBank data- 
base entries by using the BLASTn and BLASTx programs 
(www.ncbi.nim.nih.gov/blast/Blast.cgi). Primer sequences 
are shown in Table 1. The identity of bat species yielding 
specific viral products was determined by amplifying and 
sequencing a fragment of the cytochrome B gene (/0). All 
visual classifications of species were confirmed except for 
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that of T. perforatus bats. There is no reference sequence 
for T. perforatus bats in GenBank. However, because the 
closest reference sequence was from T. nudiventris bats, at 
84% identity we presume that the product represents bona 
fide T. perforatus bat cytochrome B gene sequence. Repre- 
sentative cytochrome B sequences have been uploaded to 
GenBank (accession nos. KF498635—-KF498641). 

Table 1 indicates the CoV genera identified by using 
individual primer sets. As anticipated, pan-CoV assays de- 
tected a- and B-CoVs. One assay specific for MERS CoV 
(9) also detected a-CoVs. This finding reinforces the need 
for sequence confirmation of PCR products. Table 2 indi- 
cates the CoV species identified with respect to location, 
sample type, and bat species. CoV sequences were ampli- 
fied from rectal swab samples or fecal pellets and from 
roost feces but not from serum, throat swab samples, or 
urine. a-CoV sequences were amplified more frequently 
than B-CoV sequences (223 vs. 4). Whereas a- and B-CoV 
sequences were amplified from CoVs from T. perforatus, 
E. helvum, and R. hardwickii bats, only alpha sequences 
were amplified from CoVs from P. kuhlii bat samples. 
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Table 2. CoVs detected in bats, Saudi Arabia* 


MERS CoV in Bats, Saudi Arabia 


No. samples tested (no. positive) Total no. 
Bat family, genus, No. Throat Fecal Roost samples, Total no. positive samples, n = 
species Location bats swab pellets Urine Serum feces n=1,003 227 (closely related CoVs)t 
October 2012 
Emballonuridae 
Taphozous perforatus Bisharuins 29 29(0) 25(2) 8(0) 22 10 (1) 94 1 B-Cov (1 MERS novel CoV) 
and 2 a-CoVs (1 bovine 
respiratory CoV, 1 Kenya bat 
CoV BtKY86) 
Pteropodidae 
Eidolon helvum Bisha town 25 25(0) 25(5) 13(0) 19 NA 82 1 B-Cov (1 Eidolon bat CoV- 
center HKU‘1) and 4 a-CoVs (4 Kenya 
bat CoV BtKY86) 
Rousettus Bisha town 3 3 (0) 3(0) 1(0) 2 NA 9 NA 
aegyptiacus center 
Rhinopomatidae 
Rhinopoma Naqi and 36 36(0) 35(0) 4(0) NA 15 (0) 90 NA 
hardwickii Old Nagi 
Rh. microphyllum Old Nagi 1 1 (0) 1 (0) NA NA NA 2 NA 
Vespertilionidae 
Eptesicus bottae Bisha ruins 1 1 (0) 1(0) 1(0) NA 32(0) 35 NA 
Pipistrellus kuhlii Bisha ruins 1 1 (0) 1 (0) NA NA NA 2 NA 
April 2013 
Rhinopomatidae 
Rh. hardwickii Greater NA NA NA NA NA 209 209 2 B-Covs (2 canine respiratory 
Bisha area (93) CoVs) and 91a-CoVs (5 canine 
CoVs, 2 Miniopterus bat CoVs, 
84 Chaerephon bat CoV) 
T. perforatus Bisha ruins NA NA NA NA NA 203 (0) 203 NA 
Vespertilionidae 
P. kuhlii Greater 9 9 (0) NA NA NA 263 277 126 a-CoVs (69 alphaCoV 
Unaizah (126) P.kuh-Spain, 3 canine CoVs, 
area 37 bat CoV P.pyg/Germany, 1 
human CoV NL63, 2 Rousettus 
bat CoV HKU10, 11 porcine 
epidemic diarrhea virus, 2 
Cardioderma bat CoVs, 1 
Hipposideros bat CoV HKU10) 
Greater 5 5 (0) NA NA NA NA NA NA 
Riyadh 
area 


*CoV, coronavirus; MERS, Middle East respiratory syndrome; NA, not applicable. 


+Based on BLASTn (www.ncbi.nim.nih.gov/blast/Blast.cgi). 


CoV sequences were amplified from 220 of 732 roost 
feces samples and 7 of 91 rectal swab samples or fecal pel- 
lets. A product obtained by PCR amplification of nucleic 
acid from a fecal pellet of a 7. perforatus bat captured in 
October 2012 in Bisha showed 100% nt identity to the hu- 
man B-CoV 2c EMC/2012 cloned from the index case-pa- 
tient in Bisha. A phylogenetic analysis of CoVs obtained 
in this study is shown in Figure 2. CoV sequences have 
been uploaded in GenBank (accession nos. KF493884— 
KF493888). 


Conclusions 

A wide range of CoV species are circulating among 
bats in Saudi Arabia. Although the prevalence of CoVs 
was high (~28% of fecal samples), MERS CoV was found 
in only | bat. A 3.5% MERS CoV infection rate (n = 29; 
95% CI 0-20%) in T. perforatus bats is low compared 
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with that for severe acute respiratory syndrome—like CoV 
in rhinolophid bats in China (10%—12.5%) but consistent 
with CoV prevalence among bats in Mexico (4). Further- 
more, the sensitivity for viral nucleic acid detection in 
samples collected in October 2012 was probably reduced 
because of failure in cold chain transport. Whereas 219 
(32%) of 675 of fecal pellets collected in April revealed a 
CoV sequence by PCR, only 8 (5%) of 148 of rectal swab 
samples or fecal pellets collected in October were positive 
by the same assays. We were unable to recover additional 
sequences beyond the 190-nt RdRp fragment represented 
in Figure 2 but are confident in the fidelity of the finding. 
First, although RdRp is a conserved portion of the CoV 
genome, there is no precedent for 100% identity of a bat 
sequence with a human MERS CoV sequence. Second, 
when this work began we did not have cultured MERS 
CoV, human MERS samples, or MERS CoV cDNA in the 
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71) KSA-287-Taphozous perforatus bat 
JX869059-Human betacoronavirus 2c EMC/2012 
KC164505-Betacoronavirus England 1 
KC667074-Human betacoronavirus 2c England-Qatar/2012 
KF186567-Middle East respiratory syndrome coronavirus-Al-Hasa 1 2013 
KF192507-Middle East respiratory syndrome coronavirus 
KC776174-Human betacoronavirus 2c Jordan-N3/2012 
KC243390-Bat coronavirus BtCoV/8-724/Pip pyg/ROU/2009 
KC243392-Bat coronavirus BtCoV/UKR-G17/Pip nat/UKR/2011 
KC243391-Bat coronavirus BtCoV/8-691/Pip nat/ROU/2009 
9 , KC522088-Pipistrellus bat coronavirus HKUS isolate 33R 
KC522088-Pipistrellus bat coronavirus HKUS isolate 32R 
KC522044-Tyloncyteris bat coronavirus HKU4 isolate 10R 
99 EF065505-Bat coronavirus HKU4-1 

HQ728482-Eidolon bat coronavirus/Kenya/KY 24/2006 
99 | AB514460-Human coronavirus HKU1-ORF 1b-PHLO7-1656 
DQ415911-Human coronavirus HKU1-N15 genotype B 
JN129834-Human coronavirus OC43-HK04-01 
AB514461-Human coronavirus OC43-ORF 1b-PHLO7-1581 
FJ938065-Bovine respiratory coronavirus-AH187 
KSA-004-Rhinopoma hardwickii bat 
84 JX860640-Canine respiratory coronavirus-K37 


MERS CoV 


95 


Figure 2. Phylogenetic tree 
showing genetic relatedness 
between coronaviruses 
identified in bat samples from 
Saudi Arabia (boldface), 
MERS coronaviruses, and 
other published coronavirus 
sequences available in 
GenBank. The  maximum- 
likelihood tree of partial RNA- 
dependent RNA _ polymerase 
gene (nt position 15068-15249 
of GenBank accession no. 
JX869059) was constructed 
by using the Tamura-Nei 
model with discrete gamma 
rate differences among sites 


FJs8e6ee Bat SARS Co-Rs672/2006 

GQ153547-Bat SARS coronavirus-HKU3-12 
99 L GQ153540 Bat SARS coronavirus HKU3-5 

JX503060-Human coronavirus 229E isolate 0349 


92 AY994055. 1| Feline infectious peritonitis virus 

0Q811787-PRCV-ISU-1 

HQ012368-Feline coronavirus-UU18 
97 FJ938061-Feline coronavirus-UU3 

HM245925-Mink coronavirus strain-WD1127 

0Q249235-Bat coronavirus-HKU2-1 


EU420138-Bat coronavirus 1A-AFCD62 
100 ») KSA-002-Rhinopoma hardwickii bat 
KSA-001-Pipistrellus-kuhlii bat 


KSA-003-Rhinopoma hardwickii bat 
AF353511 PEDV strain CV777 


AY567487-Human coronavirus NL63 


laboratory at Columbia University where samples were 
removed directly from the tubes in which they were col- 
lected in the field for nucleic acid extraction, PCR, and 
sequence analysis. Third, the only MERS-positive signal 
was obtained in PCR analysis of the 7. perforatus bat cap- 
tured in Bisha near the home and workplace of the MERS 
index case-patient used to generate the human B-CoV 2c 
EMC/2012 sequence. 

Bats are reservoirs of several viruses that can cause hu- 
man disease, including rabies, Hendra, Nipah, Marburg, se- 
vere acute respiratory syndrome CoV, and probably Ebola 
viruses (//—/4). Cross-species transmission from bats to 
humans can be direct, through contact with infected bats or 
their excreta, or facilitated by intermediate hosts (/5). Bat 
CoVs are typically host specific; however, MERS-related 
CoVs have reportedly been found in many bat families, in- 
cluding Vespertillionidae, Molosidae, Nyteridae, and now 
Emballonuridae (sheath-tailed bats) in Africa, the Ameri- 
cas, Asia, and Europe. We sampled only a small sample 
of bats in Saudi Arabia. Nonetheless, given the rarity of 
MERS CoV sequences detected by our survey and the 
broad distribution of MERS cases throughout the Middle 
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86 | KC175339. 1| Canine coronavirus strain 171 complete genome 


HQ728485-Miniopterus bat coronavirus/Kenya/KY 33/2006 


HQ728486-Chaerephon bat coronavirus/Kenya/KY 22/2006 


as implemented in MEGA 5.2 
(www.megasoftware.net). Each 
branch shows the GenBank 
accession number followed by a 
brief description of the sequence 
used. Scale bar indicates 
nucleotide substitutions per site. 
MERS, Middle East respiratory 
syndrome; CoV, coronavirus; 
SARS, severe acute respiratory 
syndrome; KSA, Kingdom of 
Saudi Arabia. 


HQ728480-Cardioderma bat coronavirus/Kenya/KY 43/2006 
JQ989266-Hipposideros bat coronavirus HKU10 isolate TT3A 
EU769557-Bat coronavirus Trinidad/1FY2BA/2007 

EU420139-Bat coronavirus HKU8 strain-AFCD77 


East, we speculate that there are probably other hosts. Fu- 
ture work should investigate additional bat and other wild- 
life species and domestic animals for CoV infection and 
potential linkage to human disease. 
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